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Abstract 

We extend our approach to incorporate the proton-proton (pp) Coulomb force into the three-nucleon 
(3N) Faddeev calculations from elastic proton-deuteron (pd) scattering to the breakup process. The main 
new ingredient is a 3-dimensional screened pp Coulomb t-matrix obtained by a numerical solution of the 3- 
^ . dimensional Lippmann-Schwinger (LS) equation. We demonstrate numerically that the proton-deuteron (pd) 

' breakup observables can be determined from the resulting on-shell 3N amplitudes increasing the screening 

radius. However, contrary to the pd elastic scattering, the screening limit exists only after renormalisation 
^ , of the pp t-matrices. 
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I. INTRODUCTION 



The long-range nature of the Coulomb force prevents the application of the standard techniques 
developed for short-range interactions in the analysis of nuclear reactions involving two protons. 
One proposal to avoid the difficulties including the Coulomb force is to use a screened Coulomb 
interaction and to reach the pure Coulomb limit through application of a renormalisation proce- 
dure iii,y,y,B. 

For elastic pd scattering first calculations, with modern nuclear forces and the Coulomb force 
included, have been achieved in a variational hyper spherical harmonic approach ^. Only recently 
the inclusion of the Coulomb force was undertaken also for the pd breakup reaction 01. In ml, 
contrary to [5|] where the exact Coulomb force in coordinate representation has been used directly, 
a screened pp Coulomb force has been applied in momentum space and in a partial wave basis. In 
order to get the final predictions which can be compared to the data, the limit to the unscreened 
situation has been performed numerically applying a renormalization to the resulting 3N on-shell 
amplitudes P, Q]. This allowed for the first time to analyze high-precision pd breakup data and 
provided a significant improvement of data description in cases where the Coulomb force plays an 
important role Q]. 

Recently we developed a novel approach to include the pp Coulomb force into the momentum 
space 3N Faddeev calculations [9] . It is based on a standard formulation for short range forces and 
relies on the screening of the long-range Coulomb interaction. In order to avoid all uncertainties 
connected with the application of the partial wave expansion, inadequate when working with long- 
range forces, we used directly the 3-dimensional pp screened Coulomb t-matrix. We demonstrated 
the feasibility of that approach in case of elastic pd scattering using a simple dynamical model for 
the nuclear part of the interaction. It turned out that the screening limit exists without the need 
of renormalisation not only for pd elastic scattering observables but for the elastic pd amplitude 
itself. 

In spite of the substantial progress in the pd breakup treatment achieved in P, Q] some impor- 
tant questions remained unanswered. One concerns directly the results of these calculations for 
two kinematically complete breakup geometries: the pp quasi-free-scattering (QFS) configuration, 
in which the not detected neutron is at rest in the laboratory system, and the space-star (SST) 
geometry, in which all 3 outgoing nucleons have the same momenta (magnitudes) in the plane 
which in the 3N cm. system is perpendicular to the incoming nucleon momentum. The theoretical 
predictions based on nuclear forces only show, that the cross sections for QFS and SST are quite 
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stable against changes of the underlying interactions, including also three- nucleon forces jlO|. At 
energies below w 20 MeV theory overestimates the SST pd cross sections by w 10%, and over- 
estimates the pp QFS cross sections by ^ 20%, respectively 0, [l^. With increasing energy the 



theoretical cross sections come close to the data 



lOj], which indicates that the pp Coulomb force is 



very probably responsible for these low energy discrepancies. However, the Coulomb force effects 
found in |f| are practically negligible for the pd SST configuration and only slightly improve the 
description of the pp QFS data [ll|, 



e tor tne pd bb 

i3,Q,Q,Q. 



This inability to understand the pp QFS and pd SST cross sections motivated us to reconsider 
the inclusion of the Coulomb force into momentum space Faddeev calculations. One main concern 
in such type of calculations is the application of a partial wave decomposition to the long-ranged 
Coulomb force. Even when screening is applied it seems reasonable to treat from the beginning 
the screened pp Coulomb t-matrix without partial wave decomposition because the required limit 
of vanishing screening leads necessarily to a drastic increase of the number of partial wave states 
involved. As an example we provide numbers for the exponential screening e~^Ti^ . Taking the 
screening radius i? = 20 fm and n = 4 requires all / < Imax = 10 partial wave states to reproduce 
the 3-dimensional pp screened Coulomb t-matrix at E^^^'' = 13 MeV. Increasing the screening radius 
to i? = 120 fm requires Imax ~ 50 which is a big numerical challenge. Even more, that would lead to 
an explosion of the number of 3N partial waves required for convergence. Another problem concerns 
the treatment of the pp Coulomb interaction in its proper coordinate, which is the relative proton- 
proton distance, throughout the calculations. Any deviation from this restriction can cause effects 
which are difficult to estimate. Both problems were accounted for in our formulation presented 
in 0] and applied there to pd elastic scattering. 

In the present paper we extend that approach to the pd breakup. Again we apply directly the 
3-dimensional screened pp Coulomb t-matrix without relying on a partial wave decomposition. In 

we demonstrated that the physical pd elastic scattering amplitude can be obtained from the 
off-shell solutions of the Faddeev equation and has a well defined screening limit. In contrast to 
elastic scattering, where the amplitude itself does not require renormalisation, in case of the pd 
breakup the on-shell solutions of the Faddeev equation are required. They demand renormalisation 
in the screening limit which can be achieved through renormalisation of the pp t-matrices. 

In section In] for the convenience of the reader we shortly present the main points of the formalism 
outlined in details in The screening limit for pd breakup is discussed in section Hill and the 
results shown in section IIV[ The summary is given in section |Vl 
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II. FADDEEV EQUATIONS WITH SCREENED PP COULOMB FORCE 

We use the Faddeev equation in the form [2] 

T|$ >= tP\^ > +tPGoT\<l> > (1) 

where P is defined in terms of transposition operators, P = P12P23 + ^'i3-P23i Go is the free 3N 
propagator, |$ > the initial state composed of a deuteron state and a momentum eigenstate of the 



proton. Knowing T|<I> > the breakup as we 
by quadratures in the standard manner 



IOC- 



1 as the elastic pd scattering amplitudes can be gained 



We use our standard momentum space partial wave basis \pqa > 

\pqa >== \pq{ls)j{X^)I{jI)J{t^)T > (2) 

and distinguish between the partial wave states \pqa > with total 2N angular momentum j below 
some value jmax- j < jmax^ in which the nuclear, Vn, as well as the pp screened Coulomb interac- 
tion, Vj^ (in isospin t = 1 states only), are acting, and the states \pqP > with j > jmax, for which 
only is acting in the pp subsystem. The states \pqa > and \pq(3 > form a complete system of 
states 

/ p'^dpq'^dq^lpqa) {pqa\ = / p^dpq^dqi^ \pqa) {pqa\ + ^ \pql3) (j)qf3\) = I . (3) 
a a p 

Projecting Eq.([T]) for r|<I> > on the \pqa > and \pq(5 > states one gets the following system of 

coupled integral equations 

{pqa\T\<^) = {pqa\tf,+,P\<^) 

+ {pqc^\4+cPGoYl j P^dp'q'^dq'\p'q'a'){p'q'a'\T\<^) 

a' 

+ {Pqa\ E / P^dp'q'^dq' \p'q'(3') {p'q'(3'\T |<I>) (4) 

13' 

{pq(3\T\^) = {pqP\t^P\<^>) 

+ {pqp\tfPGoJ2 / P'^dp'q'^dq'\p'q'a'){p'q'a'\T\^) 

a' 

+ {pqP\ t^PGo I P'^dp'q'^dq' \p'q'f3') {p' q' (5' \T \^) (5) 
/3' '' 

where t^^^ and are t-matrices generated by the interactions Vat + and V^, respectively. 
Namely for states |q > with two-nucleon subsystem isospin t = 1 the corresponding t-matrix 
element < pa\t^j^^{E — -^q^)\p'a' > is a linear combination of the pp, t^+c, and the neutron- 
proton (np), tnpi t = 1 t-matrices, which are generated by the interactions Vpp"'^^ + and 



y strong ^ respectively. The coefficients of tliat combination depend on the total isospin T and T' of 



states la > and la' > 
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<t = oT = ^|4+jt' = or' = ^> = (7) 

In case of only the screened pp Coulomb force Vj^ is acting. 

The third term on the right hand side of ([5]) is proportional to < pqP\t^PGo\p'q' (3' >< p'q'(3'\t^. 
A direct calculation of its isospin part shows that independently from the value of the total isospin 
T it vanishes [9]. 

Inserting < pq(3\T\^ > from ^ into ^ one gets 

{pqa\ T 1$) = {pqa\ 4^,P |$) + {pqa\ 4+,PGotf P |1>) 

- {pqa\4^,PGoY, J P'^dp'q'^dq' \p'q'a') {p'q'a'\tfP\^) 

a' 

+ {pqa\t§+^PGoY, j p'^dp'q'^dq'\p'q'a') {p'q'a'\T\^) 

a' 

+ {pqa\4^,PGot^PGoY^ j P'^dp'q'^dq' \p'q'a') {p'q'a'\Tm 

a' 

- {pqa\4^,PGoY, j p'^dp'q'^dq' \p'q'a') {p'q'a'\t^PGo 

a' 

I p"^dp"q"^dq"\p"q"a"){p"q"a"\T\<^) . (8) 

This is a coupled set of integral equations in the space of the states |a > only, which incorporates 
the contributions of the pp Coulomb interaction from all partial wave states up to infinity. It can 
be solved by iteration and Pade summation [9, llOll . 

□ 

When compared to our standard treatment without screened Coulomb force [10|] there are two 
new leading terms < pqa\t^j^^PGQt^P\(^ > and -< pqa\t§_^^PGo\a' >< a'\t^P\^ >. The first 
term must be calculated using directly the 3-dimensional screened Coulomb t-matrix t^, while 
the second term requires partial wave projected screened Coulomb t-matrix elements in the |a > 
channels only. The kernel also contains two new terms. The term < pqa\t^_^^PGot^ PGo\a' >< 



a'|T|$ > must again be calculated with a 3-dimensional screened Coulomb t-matrix while the 
second one, -< pqa\t^j^^PGQ\a' >< a'\t^PGo\a" >< a"|r|<I> >, involves only the partial wave 
projected screened Coulomb t-matrix elements in the \a > channels. The calculation of the new 
terms with the partial wave projected Coulomb t-matrices follows our standard procedure. Namely 
the two sub kernels t^j^^PGo and t^PGo are applied consecutively on the corresponding state. The 
detailed expressions how to calculate the new terms with the 3-dimensional screened Coulomb t- 
matrix are given in Appendix A of Ref. . 

The transition amplitude for breakup < $o I I ^ > is given in terms of T | by [lOl . Ill 

{^o\Uo\^) = {^o\{l + P)T\^) (9) 

where |<5o >= \pqTnim2'm3i'ih'2i^3 > is the free state. The permutations acting in momentum-, 
spin-, and isospin-spaces can be applied to the bra-state < (j)o\ =< pqmim2m3i'ii'2i^3\ changing 
the sequence of nucleons spin and isospin magnetic quantum numbers nii and and leading to 
well known linear combinations of the Jacobi momenta p, q. Thus evaluating ([9]) it is sufficient 
to regard the general amplitudes < pqmim2m3iyiiy2J^3\T \^) = {pq\ T\^). Using Eq. ([T|) and the 
completness relation ^ one gets: 



T 1$) = {pq\Y^ j p''^dp'q''^dq' \p'q'a') {p'q'a'\T\^) 



1 p'^dp'c 


1 dq 


\pqa) 


/ 1 1 i\ 
{pqa\ 


1 p'^dp'c 


fdq' 


\p'q'a') 


{p'q'a'\ 



+ {PQ I tcP\^) + {pq I t^PGo j p'^dp'q'^dq' \p'q'a') {pq'a'\T\<^) 



(10) 



It follows, that in addition to the amplitudes < pqa\T\^ > also the partial wave projected 
amplitudes < pqa\t^P\^ > and < pqa\t^PGo\a' >< a'|T|$ > are required. The expressions for 
the contributions of these three terms to the transition amplitude for the breakup reaction are 

Q 

given in Appendix B of Ref. [9|] . 

The last two terms in (jlOh again must be calculated using directly 3-dimensional screened 
Coulomb t-matrices. In Appendix C of Ref. {ol the expression for {pq\ t^P\^) and in Appendix 
D of Ref. S the last matrix element < pq \t^PGo\a' >< a'|T|<I> > are given. 
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III. THE SCREENING LIMIT 



The set of coupled Faddeev equations dSD is well defined for any finite screening radius. It is an 
exact set assuming that the strong NN t-matrix can be neglected beyond a certain jmax, which is 
justified. Further the pp screened Coulomb force is taken into account to infinite order in the partial 
wave decomposition in form of the 3-dimensional screened Coulomb t-matrix tpp. The important 
challenge is to control the screening limit for the physical pd breakup amplitude Q. In case of 
elastic scattering we provided in Ref. [9] analytical arguments and showed numerically that the 
physical elastic pd scattering amplitude itself has a well defined screening limit and does not require 
renormalisation. This can be traced back to the fact that to get the elastic pd scattering amplitude 
it is sufficient to solve the Faddeev equations dH]) for off-shell values of the Jacobi momenta 

^ + ±q^^E. (11) 
m 4m 

The off-shell Faddeev amplitudes {pqa\ T\^) of Eq.® are determined by off-shell nucleon-nucleon 
t-matrix elements t{p,p';E — 4^9^), which have a well defined screening limit (see the following 
discussion and examples). 

Contrary to pd elastic scattering the physical breakup amplitude ([9]) corresponds to the on-shell 
values of Jacobi momenta 

^ + ^q' = E^ ^ql^^ . (12) 
m 4m 4m 

That means that the physical pd breakup amplitude ([TO]) requires on-shell Faddeev amplitudes 



{poqa\ T |<I>) together with the four, also on-shell, additional terms in dSl), with po = \J\{qmax ~ q'^^ 
The on-shell Faddeev amplitudes can be obtained from the off-shell solutions {pqa\ T |<I>) using 

{poqa\ T 1$) = {poqa\ 4+^^ |«>) + {poqa\ 4^,PGotfP \<P) 

- {Poqa\ 4+cPGo j P'^dp'q'^dq' \p'q'a') {p' q' a' P \<l>) 

a' 

+ {poqa\4+,PGoJ2 J p'^dp'q'^dq' \p'q'a') {p'q'a'\T \<^) 

a' 

+ {poqa\4+,PGotfPGoY^ f p'^dp'q'^dq' \p'q'a') {p'q'a'\T\^) 



- {poqa\4_,,PGoY, J p'^dp'q'^dq' \p'q'a') {p'q'a'\t^PGo 

a' 

[ p"^dp"q"^dq"\p"q"a"){p"q"a"\T\<^) . (13) 

These on-shell amplitudes together with additional, also on-shell, terms in (jlOh define the physical 

2 

breakup amplitude ([9]). That in consequence requires half-shell t-matrix elements t{po,p'; ^) which 



are of 3 types: the partial wave projected pure screened Coulomb generated by V^, the partial 
wave projected generated by Vstrong + Vj^, and the 3-dimensional screened Coulomb matrix 
elements. 

It is well known Q, [3] that in the screening limit R ^ oo such half-shell t-matrices acquire 
an infinitely oscillating phase factor e**^^^\ where ^r{p) depends on the type of the screening. 
For the exponential screening its form depends on two parameters, the screening radius R and the 
power n: 

= ^e-(*)" . (14) 
At a given value n the pure Coulomb potential results for R oo. As has been shown in [2^ 



based on 
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22l |. the related phase ^r{p) is given as 

= -r][ln{2pR) - e/n] (15) 

where e = 0.5772 ... is the Euler number and rj = the Sommerfeld parameter. 

Contrary to the half-shell, the off-shell t-matrix elements do not acquire such an oscillating 
phase and their screening limit is well defined. 

In Figs. [33] we demonstrate that behavior for the 3-dimensional half-shell screened Coulomb 
pp t-matrix [23] . Increasing the screening radius R changes drastically the imaginary part of the 
t-matrix (Fig. [2^). The real part is more stable but does not approach the pure Coulomb limit 
(Fig. [1^). Renormalizing by the phase factor e"**-^^^^ of Eq. (I15|) provides a well defined limit to 



the pure Coulomb half-shell result of Ref. [24f] (Fig. [TJa and [2)3). 

For the 3-dimensional off-shell screened Coulomb pp t-matrix the pure Coulomb screening limit 
of Ref. [25, 



26( 1 is achieved without any renormalisation factor for screening radia i? > 20 fm 



(Fig. [3D. 

Analogous behavior for the partial wave decomposed / = half-shell screened Coulomb and 
the ^5*0 t^+c t-matrices is shown in Figs. [D and [5[ respectively. While the imaginary part again 
exhibits drastic changes when the screening radius increases (Fig. [4^ and[5K), removing the phase 
factor e~**^*^'P) (renormalisation) provides a well defined limit for the screening radia i? > 40 fm 
(Fig. [3)3 and[S)3). It is seen that in case when the screened Coulomb potential is combined with the 
strong force also the real part of the half-shell t-matrix undergoes strong changes with increased 
screening (Fig. [5^). 

For the partial wave decomposed off-shell I = screened Coulomb and the ^Sq t-matrix 
elements a well defined screening limit is reached without any renormalisation (Fig. [6j3 and [6^, 
respectively) . 



That oscillatory phase factor appearing in the half-shell proton-proton t-matrices requires a 
carefull treatment of (|13p to get the screening limit for the {poqa\T\^) amplitudes. Namely for 
the states \a > with the two-nucleon subsystem isospin t = 1 the corresponding t-matrix element 

2 

< PoC(\tN+c(^)\p''^' > is a linear combination of the pp and neutron-proton (np) t = 1 t-matrices, 
the coefficients of which depend on the total isospin T and T' of the states \a > and \a' > (see 
discussion after dSJ). It follows that to achieve the screening limit one needs to renormalize the 
pp t-matrix tpp_^_^ in that combination before performing the action of the operators in (jl3p . The 
term in that linear combination coming with the np t-matrix tnp does not require renormalisation. 



IV. NUMERICAL RESULTS 



To demonstrate the feasibility of our approach we applied the outlined formalism to a simple 
dynamical model in which the nucleon-nucleon force was restricted to act in ^Sq and '^S'l — Di 
partial waves only and taken from the CD Bonn potential j^]- The proton-proton Coulomb force 
was modified by the exponential screening (jl4p with the screening radius R and n = 1. 

To investigate the screening limit R ^ oo we generated set of partial-wave decomposed t- 
matrices, t^, based on the screened pp Coulomb force only or combined with the strong pp inter- 
action, tpp_^_^, taking R = 20, 40, 60, 80, 100, 120 and 140 fm. With that dynamical input we solved 
the set of Faddeev equations ([8]) for off-shell values of the Jacobi momenta p and q and for the total 
angular momenta of the p-p-n system up to J < ^ and both parities. Then the on-shell Faddeev 
amplitudes {poqa\ T |$) were gained through ([13]). In this first study we restricted ourselves to the 
perturbative approximation for the 3-dimensional screened Coulomb t-matrix: = V^. Of course 
in the future studies that approximation will be avoided and the full solution of the 3-dimensional 
LS equation for the screened pp Coulomb t-matrix will be used j23i] . When calculating observ- 
ables we also omitted the last term (^>o|(l + P)t^PGoT\^^ in pUj) coming with the 3-dimensional 
screened Coulomb t-matrix. 

The results for the breakup reaction are shown in Figs. [TlfTOl where the exclusive cross sections 
dQi^hdS QFS and SST configurations parametrized through the arc-length of the kinematical 

S-curve are presented. 

For the QFS and SST (see Fig. [7] and [H respectively) the convergence in the screening radius is 
achieved at i? = 60 fm. For QFS the Coulomb force decreases the cross section with respect to the 
nd case and brings the theory close to the pd data. For SST the Coulomb force also brings theory 
close to the pd data, however, only at S-values close to the space-star condition {S ^ Q MeV). For 
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S-values further away the theory is far above the pd data. 

The theoretical prediction for both geometries results through interference of different terms 
contributing to the breakup amplitude. The importance and magnitudes of the contributions 
coming from different terms in the breakup amplitude differs for those two geometries (see 
Fig. [5] and Fig. [TU|) . In both cases the largest is the contribution of the first term in pO|) 
($o|(l + P)\a >< a\T\^) (black dashed-dotted line in Figs. [8] and [JO]). For QFS and SST the cross 
section resulting from that term is below the pd data and below the full result which encompasses all 
terms (solid line). The magnitudes of three additional terms: (^'oKl + P)t^P\^) (the fourth term 
in (jlOp calculated with the 3-dimensional screened Coulomb t-matrix, here approximated by Vj^, 
and given by the green short-dashed line), the second term in (jlOp (<I>o|(l + >< a|t^P|$) 
(calculated with the partial- wave projected screened Coulomb t-matrix and given by the blue 
dashed-double-dotted line), and the third term in ([10]) (<I>o|(l + P)|a >< a\t^PGoT\^) (calcu- 
lated again with the partial- wave projected screened Coulomb t-matrix and given by the maroon 
double-dashed-dotted line), are small. Because they are difficult to see on the scale of Figs. [8^ and 
Fig. [TOk they are again presented in the part b) of these figures. For both configurations the term 
with the 3-dimensional t-matrix (<I>o|(l + P)i!^P\^) gives the smallest contribution. Smallness of 
these terms does not mean however, that they are unimportant because the interference effects are 
nonnegligible and act in different ways for QFS and SST. 



double-dashed-dotted line). For QFS and SST taking the amplitude of that second largest contri- 
bution together with (<I>o|(l + P)\oi >< a|r|<I>) changes significantly the cross section. For QFS it 
is the black dotted line in Fig. [HJi resulting from (^>o|(l + P)\a >< a\{T - t^P)\<^) while for SST 
it is blue long-dashed line in Fig. [TUh- resulting from (<I>o|(l + P)\a X a\{T - t^PGoT)\^). 

The third largest contribution provides for both configurations smaller changes of the 
cross section and the result when the second and third largest contributions are included 
($o|(l + P)\a X a\{T - t^PGoT - t^P)\^) is given by the red long-dashed line. It is above 
the pd data for both geometries. 

Finally, including the smallest contribution from the 3-dimensional screened Coulomb t-matrix 
(<I>o|(l + P)t^P\^) brings the theory to the pd data for the QFS geometry at all S-values and for 
the SST configuration at S-values close to the space-star condition. 
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V. SUMMARY AND CONCLUSIONS 



We extended our approach to include the pp Coulomb force into the momentum space 3N 
Faddeev calculations presented in Ref. (ol for elastic scattering to the pd breakup. It is based on a 
standard formulation for short range forces and relies on the screening of the long-range Coulomb 
interaction. In order to avoid all uncertainties connected with the application of the partial wave 
expansion, unsuitable when working with long-range forces, we apply directly the 3-dimensional 
pp screened Coulomb t-matrix. 

Using a simple dynamical model for the nuclear part of the interaction we demonstrated fea- 
sibility of that approach for the treatment of the pd breakup. We demonstrated that contrary 
to the pd elastic scattering, where the resulting amplitudes do not require renormalisation, it is 
unavoidable to perform renormalisation of the pp half-shell t-matrices in order to get the physical 
breakup amplitude. Namely that amplitude has two contributions, one driven by the interaction 
in the pp subsystem and second in the np subsystem. Only the first part requires renormalisation. 

We have shown that converged results for breakup can be achieved with finite screening radia. 

We calculated contributions of different terms to the breakup cross section in QFS and SST 
configurations. The interference between different contributions leads to an interference pattern, 
which is different for QFS and SST configurations. In our restricted dynamical model the pp 
Coulomb interaction brings the nd breakup cross sections close to the pd data for the QFS con- 
figuration. Also for the SST geometry in the vicinity of the space-star condition the pd theory is 
close to the pd data. However, further away on the S-curve the theory lies above the data. 

In future studies the perturbative approximation for the 3-dimensional screened Coulomb t- 
matrix will be avoided and higher partial wave components of the nucleon-nucleon interaction will 
be included. 
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FIG. 1: (color online) The real part of the nonrenormalized (a) and renormalized (b) 3-dimensional half- 
shell screened Coulomb pp t-matrix tf{po,p', x;^). The lines correspond to the exponential screening with 
n = 1 and different screening radia: i? = 20 fm (black dotted line), i? = 60 fm (blue dashed-dotted line), 
i? = 120 fm (red dashed Hue), i? = 180 fm (black solid line). The Coulomb half-sheU result of Ref. [2J] is 
given by thick dots. The momentum pq — 0.396 fm^^ and x ~ 0.706. 
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FIG. 2: (color online) The imaginary part of the nonrenormalized (a) and renormalized (b) 3-dimensional 
half-shell screened Coulomb pp t-matrix tf {pQ,p' ,x; ^). For the description of the lines and values of po 
and X see Fig. [TJ 
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FIG. 3: (color online) The 3-dimensional off-shell screened Coulomb pp t-matrix tf{p,p',x] — The lines 
correspond to the exponential screening with n = \ and different screening radia: i? = 20 fm (black dotted 
line), i? = 60 fm (blue dashed-dotted line), R — 120 fm (red dashed line), R — 180 fm (black solid line). 



The Coulomb off-shell result of Ref. 
p = 0.375 fm-i and x = 0.706. 



25 



26[ is given by thick dots. The momentum po = 0.396 fm 
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FIG. 4: (color online) The real (upper panels) and imaginary (lower panels) parts of the nonrenormalized 

2 

(a) and rcnormalizcd (b) I = half-shell screened Coulomb pp t-matrix t^{po,p', ^). The lines correspond 
to the exponential screening with n = I and different screening radia: i? = 20 fm (black dotted line), 
i? = 40 fm (blue short-dashed line), R — 60 fm (brown long-dashed line), i? = 80 fm (red short-dashed- 
dotted line), R = 100 fm (maroon long-dashed-dotted line), R = 120 fm (green short-dashed-double-dotted 
line), R — 140 fm (blue solid line). The momentum po — 0.26 fm^^. 
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FIG. 5: (color online) The real (upper panels) and imaginary (lower panels) parts of the nonrenormalized 

2 

(a) and renormalized (b) ^5*0 half-shell pp t-matrix tpp^^{pQ,p' , ^). The lines correspond to the exponential 
screening with n = 1 and different screening radia: i? = 20 fm (black dotted line), i? = 40 fni (blue short- 
dashed line), i? = 60 fm (brown long-dashed line), i? = 80 fm (red short-dashed-dotted line), i? = 100 fm 
(maroon long-dashed-dotted line), R ~ 120 fm (green short-dashed-double-dotted line), R — 140 fm (blue 
solid line). The momentum po = 0.26 fm^^. 
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FIG. 6: (color online) The real (upper panels) and imaginary (lower panels) parts of the tpp+ciP^P' ^ m) 

2 

(a) and the I = screened Coulomb t^{p,p', ^) (b) ofF-shell t-matrices. For the description of the lines see 
Fig. m The momentum po — 0.26 fm^^ and p — 2.38 fm^^. 
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FIG. 7: (color online) The convergence in the cut-off radius R of the d{p,pip2)n breakup cross section in 
a kinematically complete QFS configuration with polar angles of the two outgoing protons 9i = 62 = 39° 
and azimuthal angle (j>i2 = 180°. The incoming proton energy is E^"''' = 13 MeV and theoretical predictions 
are based on a screened Coulomb force and the CD Bonn nucleon-nucleon potential [27I restricted to ^5*0 
and '^Si-^Di partial waves. The screening radius is: i? = 20 fm (black dotted line), i? = 40 fm (green 
dashed-double-dotted line), i? = 60 fm (blue dashcd-dotted line), i? = 80 fm (red double-dashed-dotted 
line), R — 100 fm (blue dashed line), R — 120 fm (red long-dashed line), R — 140 fm (black solid line). The 
black long-dashed-dottcd line is the break-up cross section with the pp Coulomb interaction switchcd-off. 
The pluses are EL'''' = 13 MeV pd breakup data of Ref. 
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FIG. 8: (color online) In part a): the contribution of different terms to the cross section of the 
QFS configuration of Fig. [7l The (black) dashed-dotted line is the contribution of the first term 
($o|(l + -P)|a >< Q;|r|$) in (jTU]) and the (green) short-dashed line is the contribution of the fourth term 

coming with the 3-dimensional screened Coulomb t-matrix (in the present 
calculation = Vj^)- The (blue) dashed-double-dotted and (maroon) double-dashed-dotted lines are 
contributions of the second ($o|(l + P)\a >< "l^f ^^1*) and of the third ($o|(l + P)\a >< a\t^PGoT\^) 
term in (jlOp . respectively, which are calculated with partial- wave decomposed screened Coulomb t-matrix. 
The (black) dotted and (blue) long-dashed lines result from the ($o|(l + ^')|a >< "1(7" - P)|$) and 
{^o\{^ + P)\oi X a\{T - t^PGoT)\^) amplitudes, respectively. The (red) long-dashed line is the contri- 
bution of the ($o|(l + P)\a >< a\(T - <f P - t^PGoT)\^) amplitude. The (black) solid line encompasses 
all four terms. All results are for screening radius R = 100 fm. The part b) of the figure shows contributions 
of small terms which are difficult to see on the scale of part a). 
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FIG. 9: (color online) The convergence in the cut-ofF radius R of the d{p,piP2)n breakup cross section in a 
kincmatically complete SST configuration with polar angles of the two outgoing protons 9i = 62 = 50.5° and 
azimuthal angle 0i2 = 120°. The incoming proton energy is E'^"'' = 13 MeV and the theoretical predictions 
are based on the screened Coulomb force and the CD Bonn nucleon-nucleon potential 
and '^Si-^Di partial waves. For the description of the lines see Fig. [71 The x-es are 
breakup data of Ref. [ill ^^'^ pluses are i?'"'' = 13 MeV nd breakup data of Ref. 28]. 
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FIG. 10: (color online) In part a): the contribution of different terms to the cross section of the SST 
configuration of Fig. [9] The solid line encompasses all terms. For the description of other lines see Fig. [8l 
All results are for screening radius R — 100 fm. The part b) of the figure shows contributions of small terms 
which are difficult to see on the scale of part a) . 
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